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We extend the quasiclassical formalism for diffusive superconductors by deriving anisotropic (gra¬ 
dient) corrections to the Usadel equation. We demonstrate that in a number of physical situations 
such corrections may play a crucial role being responsible for the effects which cannot be recovered 
within the standard Usadel approximation. One of them is the so-called photoelectric effect in 
superconductors and superconducting-normal (SN) hybrid structures. Provided a superconducting 
part of the system is irradiated by an external ac electromagnetic field the charge imbalance devel¬ 
ops and a non-vanishing dc voltage is induced across the SN interface. In the presence of magnetic 
impurities in a superconductor the magnitude of this effect becomes large and can easily be detected 
in modern experiments. 
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I. INTRODUCTION 

Quasiclassical theory of superconductivityic— is a pow¬ 
erful tool for comprehensive analysis of the supercon¬ 
ducting state. Quasiclassical methods are well suited 
to study a broad range of phenomena with character¬ 
istic energies/frequencies well below the Fermi energy ep 
and characteristic lengths strongly exceeding the Fermi 
wavelength or, equivalently, the inverse Fermi momen¬ 
tum 1/pf- Within this approach all relevant information 
about the superconducting state can be extracted from 
the energy-integrated matrix Green-Eilenberger function 
g which generally depends on the electron momentum di¬ 
rection as well as on both coordinate and time variables. 

The quasiclassical Eilenberger formalism also accounts 
for electron scattering at various impurities and bound¬ 
ary imperfections which are inevitably present in virtu¬ 
ally any metallic sample. In general such scattering pro¬ 
cesses can make the quasiclassical electron trajectory in 
a metal very complicated in which case the function g 
may become rather involved. Usadel^ demonstrated that 
substantial simplification can be achieved provided the 
concentration of non-magnetic impurities in a system is 
sufficiently large. In this case electron scattering at such 
impurities turns its trajectory into a diffusive one and the 
matrix function g becomes almost isotropic in the mo¬ 
mentum space. The Eilenberger equation then reduces 
to the Usadel one for the momentum averaged quasiclas¬ 
sical Green function matrix^ii^. 

Usually the Usadel approximation works well in the so- 
called dirty limit, i.e. provided the electron elastic mean 
free path i remains much shorter than the superconduct¬ 
ing coherence length of a clean system fo- However, there 
exists a number of interesting phenomena which cannot 
be captured within the Usadel approximation even in 
the limit £ <C Co- One of them is the thermoelectric ef¬ 
fect in superconductors doped with magnetic impurities^ 
as well as in superconducting hybrids with spin-active 
interfaces'^— . Different scattering rates for electrons and 


holes at magnetic impurities or interfaces yield electron- 
hole imbalance generation in such systems which in turn 
results in a large enhancement of thermoelectric cur¬ 
rents. This electron-hole imbalance effectively implies 
anisotropy of the quasiparticle distribution function in 
the momentum space^^— . It turns out that this effect 
cannot be recovered from the standard Usadel equation 
for the momentum-averaged quasiclassical Green func¬ 
tions. One of the aims of this paper is to extend further 
the Usadel perturbative expansion in the parameter 
and to derive the leading order correction to the Usadel 
equation which would enable one to properly account for 
the momentum anisotropy of the distribution function 
mentioned above. 

Another our aim is to employ this formalism in 
order to study the so-called photoelectric effect in 
superconducting-normal (SN) hybrid structures. The 
essence of this effect is the appearance of a dc voltage 
drop across a superconductor as a result of its inhomo¬ 
geneous absorption of light. It is commonly accepted 
that an external electromagnetic ac field cannot cause 
any charge imbalance in a superconductor—: Photon ab¬ 
sorption just produces a pair of electron-like and hole-like 
excitations and does not lead to any shift of the chemical 
potential. This observation remains applicable as long 
as both the density of states and quasiparticle scattering 
rates are symmetric with respect to electrons and holes. 
Electron-hole imbalance occurs provided one takes into 
account the energy dependence of the normal state den¬ 
sity of states near the Eermi energyiS. The correspond¬ 
ing effect is relatively small for generic metals without 
anomalies in the density of states in the vicinity of the 
Fermi level. Accordingly, the magnitude of the photoelec¬ 
tric effect in conventional superconductors is expected to 
be small in the parameter Tdep, where is a critical 
temperature of a superconductor. 

The situation changes dramatically if one considers 
superconductors doped with randomly distributed mag¬ 
netic impurities. In this case scattering of quasiparticles 
on such impurities may lead to significant electron-hole 
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asymmetry resulting, in turn, in a strong enhancement of 
the photoelectric effec t^^’^^ which amplitude is not any¬ 
more controlled by the small parameter Tdep- Interest¬ 
ingly enough, it turns out that this large photoelectric 
effect cannot be captured within the usual Usadel equa¬ 
tion formalism. In order to correctly describe this effect 
in diffusive superconductors doped with magnetic impu¬ 
rities it is necessary to go beyond the standard Usadel 
approximation and to inclnde anisotropic (gradient) cor¬ 
rections. This will be demonstrated below in our paper. 

The structure of the paper is as follows. In section 
El we will extend the Usadel expansion of the quasiclas- 
sical Green-Eilenberger functions and derive the leading 
order correction to the Usadel equation describing super¬ 
conductors in the diffusive limit. In section IIIII we will 
demonstrate the advantages of our approach by briefly re¬ 
deriving the expression for large thermoelectric currents 
which may flow in superconductors doped with magnetic 
impurities. In section m we will apply our formalism 
in order to analyze the photoelectric effect in SN struc¬ 
tures with magnetic impurities. Our main conclusions 
are briefly presented in section V. Some technical details 
of our calculation are displayed in Appendix. 


where (Tq is the unity matrix in the spin space and A is 
the BCS order parameter in our system. The two self¬ 
energies Em and Enm account for electron scattering on 
respectively magnetic and non-magnetic impurities con¬ 
tained in our superconductor. In the case of isotropic 
elastic impurity scattering the self-energy Enm reads 

Snm = {g ), (6) 

where angular brackets (• ■ ■) denote averaging over the 
Fermi surface. The expression for Em will be specified 
below. For the sake of simplicity we will assume that the 
matrix Q is independent on the momentum direction pp. 

Electric current density j in our system is determined 
via the Keldysh part of the quasiclassical Green function 

j = j de {vf Sp [g^h]) , (7) 

where e is the electron charge and Nq is the normal den¬ 
sity of states at the Fermi level. 


II. QUASICLASSICAL DESCRIPTION OF 
DIFFUSIVE SUPERCONDUCTORS 


A. Eilenberger formalism 


Our starting point is the quasiclassical Eilenberger 
equation which can be expressed in the form 

ivpdrg -k [fi - Enm, ff] = 0, (1) 

where vp = pp /m is the Fermi velocity, dr denotes spa¬ 
tial gradient and g{e,t,pF,r) is the quasiclassical ma¬ 
trix propagator which has the standard structure in the 
Keldysh space 


g = 




and obeys the normalization condition 




= 1 . 


( 2 ) 

( 3 ) 


B. Expansion in spherical harmonics 

Let us now assume that the electron elastic mean free 
path in a superconductor is shorter than its “clean limit” 
coherence length f <C Co vp/ max(|A|,r). In this case 
the electron motion becomes diffusive and the Green- 
Eilenberger function g only weakly depends on the mo¬ 
mentum direction piT. In order to account for this depen¬ 
dence one can formally expand the function g in a series 
of spherical harmonics 


g — go + gi^i 


gij 


(^iUj 



( 8 ) 


where n = Pf/pf is the unity vector in the direction of 
the quasiparticle momentum and the matrices go, gi, gij 
do not depend on the electron momentum. Here and be¬ 
low the summation over repeated indices is implied. As 
it will become clear below, an effective expansion param¬ 
eter in Eq. ® is VWo, i-e. 


Here and below the “check”-symbol denotes 8x8 ma¬ 
trices in KeldyshG)Nambui8)Spin space while the“hat”- 
symbol labels 4x4 matrices in Nambu®Spin space. The 
brackets in Eq. © and below imply the commutator 
[a, &] = ab — ba. The matrix Cl has the form 

ii = efa -|- A — Em, (4) 


where e is the quasiparticle energy, fa is the Pauli ma¬ 
trix in the Nambu space and A is the order parameter 
matrix which contains only “retarded” and “advanced” 
components, 



f 0 Acro\ 
l^-AVo 0 ) ’ 



The standard Usadel approximation^ amounts to keep¬ 
ing only the first two terms in the expansion ([8]). Under 
this approximation the Eilenberger equation © reduces 
to the Usadel one^i^. Here we go beyond this approx¬ 
imation and also keep the next order term in Eq. (|S]) 
which contains the matrix pij. Without any loss of gen¬ 
erality, we can assume that this matrix gij is symmetric 
and traceless, i.e. 


A = 


( 5 ) 


9ij — 9jii 9ii — 


( 10 ) 
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In order to proceed it is convenient to rewrite gradient 
term in the Eilenberger equation (ED in the operator form 
as drQ = Vg — g\7. Then Eq. ED reads 

[ii + iVpy - Enm, g] = 0. (11) 

Substituting the expansion © into Eq. EID and collect¬ 
ing all terms associated respectively with harmonics 1, 
rii and riirij — 6ij /3, we arrive at the following system of 
equations 


[^,9o] + ^ [^i,9i\ = 0, 




-I- ivF [Vj, 5o] + [Vfe, 5ife] = 0, 


[go.gik] +£ ( [Vi,5fc] -f [Vk,gi] - -^Sik [V„,p„]) = 0. 


( 12 ) 

(13) 

0 . 

(14) 


Analogously from the normalization condition ([3]) we ob¬ 
tain 


ffo + -^Mi = 1, 


(15) 


+ Mo + riMgk + h9ik) = 0 . ( 16 ) 
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9o9ik + 9ik90 T TT (9i9k + 9k9i “ 6. (17) 


For the sake of brevity here we neglected the vector po¬ 
tential which is formally set equal to zero A = 0. If 
needed, the vector potential can easily be restored by 
gauge invariance arguments implying the replacement 


V V +i-A{r,t)T 3 , 
c 


(18) 


where c is the speed of light. 

In order to evaluate gij it is sufficient to find gt only in 
the lowest order in the parameter £/^o- From Eqs. (fT31) 
and ESD we get 


9i = -£[9o^i9o - Vi) (19) 


C. Correction to Usadel equation 

From Eqs. ESD and EID we observe that the normal¬ 
ization condition for the isotropic part of the function g^ 
acquires a nontrivial correction 

-2-1 _i- - -1 

9o 1 ^9i9i — 1 2 

- 5oVi5oVi - Vi^oViffo)- (21) 

At this stage it is convenient for us to introduce the func¬ 
tion G related to go as 

go = G-j (GV.V, + V,V,G - 2V,GV,) . (22) 

It is easy to verify that this function obeys the standard 
normalization condition 


(^oViVi^o + ViVi- 


= 1. (23) 

Combining Eqs. ESD, ESD, ESD and (1221) we also express 
the matrix gi via G and hnd the correction to this matrix. 
The corresponding rather lengthy expression is displayed 
in Appendix, see Eq. (ED- Substituting this expres¬ 
sion into Eq. EID we arrive at the result (jA2l) which 
extends the Usadel equation^i^ by including terms de¬ 
scribing weak anisotropy effects in the momentum space. 

Though Eq. (IA2I) still looks complicated, in the ma¬ 
jority of cases it can be simplified further. For instance, 
in the absence of the magnetic field (A = 0) and pro¬ 
vided the matrix G depends only on one coordinate (x) 
Eq. (1221) reduces to 


9o = G- jdlG, 


while Eq. 


becomes 


iDd,{Gd^G) = [n,G] 




(24) 


(25) 


which is just the well known expression for the ampli¬ 
tude of the first harmonics within the standard Usadel 
approximation^. Combining this expression with Eqs. 
El and EZD, we obtain 


9ik — 
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ViVfego + VfeVigo + 5oViVfc + goVfcVi-l- 


+ VjgoVfc -I- VfcgoVj - dgoViffoVfego - ^go'^kgo'^^go+ 


where D = is diffusion coefficient. The first line of 

Eq. (l25|) represents the standard Usadel equation while 
its second line defines the leading order gradient correc¬ 
tions to this equation. 

The a;-component of the current density reads (cf. Eq. 

(ESD) 

7- = - Sp (g^Ts) , (26) 


U ^ ^ik (3^0 Vn^O V ngo Vn^O V u 


^ ngo 90^ n) ■ 

( 20 ) 


where aN = 2e^NoD is the normal state Drude conduc¬ 
tivity and g^ is the Keldysh component of the matrix 


For the sake of generality here and below we make a 
distinction between the products and VfcVi being 

different from each other in the presence of the magnetic 
held. 


g, = -tGd^G 


i— [n,d,G] 


-D 


e 

20vf 


{3[G,dlG] + 7[d,G,dlG]). 


( 27 ) 
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The expression (l27l) follows directly from Eq. 

It is important to emphasize that the above results 
remain applicable as long as significant changes of the 
function G occur at distances exceeding the electron elas¬ 
tic mean free path £. Should this condition be violated 
(as it is the case, e.g., in diffusive hybrid structures with 
spin-active interfaces^) the whole approach based on the 
expansion m may become insufficient. On the other 
hand, usually the function G in diffusive superconductors 
changes at distances of order of the ’’dirty limit” coher¬ 
ence length ^ ^ (which also follows from the Usadel 
equation) or even at longer distances. Hence, the effec¬ 
tive expansion parameter in Eq. ® is as 

we already indicated above in Eq. ®. Further simplifi¬ 
cations occur provided G changes at very long distances, 
in which case the terms with higher order gradients in 
Eqs. (ESI) and (l27)) can be neglected. 


III. THERMOELECTRIC CURRENT 

Let us now turn to practical calculations employing the 
above formalism. Our first example is the thermoelectric 
effect in superconductors. If a non-vanishing tempera¬ 
ture gradient VT is applied to a superconducting sample 
quasiparticle distribution function is driven out of equi¬ 
librium and the quasiparticle current develops in the sys¬ 
tem. As a result, the total current density consists of two 
contributions: 


j{r) = js{r) + a{r)drT{r), (28) 


where jsi'f') defines the supercurrent density which we 
are not interested in here and the last term represents the 
thermoelectric current. For simplicity below we restrict 
ourselves to the linear response regime, i.e. we assume 
that both the supercurrent and the thermoelectric cur¬ 
rent remain small as compared to the critical (depairing) 
current. In this case the thermocurrent is proportional 
to VT and a is the so-called thermoelectric coefficient. 

Without loss of generality we can orient the a;-axis in 
the direction of non-zero temperature gradient, i.e. we 
assume T = T{x). Employing Eqs. (051) . (E51) and (l?7l) 
we directly evaluate the thermoelectric current. Provided 
the temperature varies at a typical length scale strongly 
exceeding the superconducting coherence length ^ the 
terms with spatial gradients in Eq. (1^51) can be safely 
neglected as they could only yield small corrections to 
the Green function proportional to {dxT)^ or d^T. As 
a result, it suffices to solve Eqs. (051) within the ap¬ 
proximation assuming local equilibrium in our supercon¬ 
ductor. Then both the Green-Keldysh function and the 
self-energy take a simple form 


G^{s, x) 
x) 


G^{s, x) 
x) 


G"^{s, x) 
x) 


tanh 


e 

2T{x) ’ 


tanh 


e 

2T{x) ’ 


(29) 

(30) 


where the retarded and advanced Green functions obey 
Eq. (1^ which now reduces to 


^’"^{e,x),G^’^{e,x) = 0 , ^G^’^{e,x)^ = 1 . 


n 

Eq. (1311) can easily be solved with the result 

qR.A 


qR,A ^ ^ 


(5K.A)2 


(31) 

(32) 


where define a traceless part of the matrix 

H^A = (iR,A _ Sp(H^’'^) = 0, (33) 

H«'^ = isp(H«’-"). (34) 

For simplicity here we assume that our system is spin 
isotropic and, hence, the matrix H^A j^g^g structure 
in the spin space. Provided the superconducting order 
parameter A is purely real, the matrix Green functions 
have the following structure 

^ / /^R,A wR,A \ 

qR^a _ [ ^ _qR,aj > 


and the combination Sp(f 3 g,^) in Eq. (1261) can be evalu¬ 
ated in a straightforward manner. As we already pointed 
out, the terms with higher derivatives in Eq. E3 can 
be safely omitted provided the temperature gradient re¬ 
mains sufficiently small. With this in mind and making 
use of Eqs. (EH)-([35]) we obtain 


(36) 


^ VF^ ° T2cosh^ [e/(2T)]’ 


where iy(£) = (1/2)[G^ — G^] is the density of states. 
Substituting this combination into Eq. (l26l) , we immedi¬ 
ately arrive at the expression for the thermoelectric co¬ 
efficient 


2^ J ^^T2 cosh^ [e/(2T)]' 


(37) 


Note that an attempt to evaluate the expression Eli) 
making use of the standard Usadel equations immedi¬ 
ately yields an incorrect result a = 0. This is because in 
local equilibrium one has Sp(f3Gi9a;G)^ = 0 and, hence, 
the contribution to the electric current from the first term 
in the right-hand side of Eq. (I?7l) vanishes identically. In 
order to recover the correct expression for the thermo¬ 
electric coefficient it is necessary to go beyond the Us¬ 
adel approximation and to include gradient corrections 
evaluated in the previous section. E.g., in the case of 
superconductors doped with randomly distributed mag¬ 
netic impurities^d^ these correction terms allow to prop¬ 
erly account for the electron-hole imbalance effects. For 
Oq Eq. (1571) yields a non-zero result which - in 
the limit ^ considered here - exactly coincides with 
the expression for a derived within a somewhat different 
approach^. The corresponding calculation is already pre¬ 
sented in Ref. d, therefore it is not necessary to go into 
further details here. 
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IV. PHOTOELECTRIC EFFECT 

Let us apply the formalism developed above to the de¬ 
scription of the photoelectric effect in superconductors 
doped with magnetic impurities. Here we will concen¬ 
trate on a hybrid structure depicted in Fig. 1. This struc¬ 
ture consists of a superconductor and a normal metal sep¬ 
arated by an insulating barrier with cross section A and 
normal state resistance Rn- We will assume that the 
superconducting part of this SN structure is uniformly 
irradiated by an ac electromagnetic field and that the 
corresponding penetration depth for this field exceeds 
the thickness of the S-film. In the presence of magnetic 
impurities a combined effect of an ac field and electron 
scattering at such impurities yields charge imbalance pro¬ 
duction in a superconducto r^^d^ which, in turn, implies 
the appearance of a nonvanishing dc voltage V across 
the insulating barrier. This effect can be detected, e.g., 
by measuring the dc current I = V/Rm across the SN 
interface. The current across the SN interface with low 
transmission can be expressed in the forir>i^ 



FIG. 1: (Color online) Superconductor/normal metal tunnel 
junction irradiated at microwave frequencies. 

presence of an ac electromagnetic field which we describe 
by means of the time-dependent vector potential 

Ec 

A{t) = —cos (jjt, (42) 

UJ 

where w is the radiation frequency and E denotes the 
amplitude of the corresponding electric field component. 
Provided the intensity of the electromagnetic field is suf¬ 
ficiently small, i.e. 


I = j de{\v:^\D{pF)S^{g^ - g^hN + hNg^)) , 

(38) 

where are retarded, advanced and Keldysh Green- 

Eilenberger functions on a superconducting part of the 
interface, hN^e) = tanh(e/2T) is related to the equilib¬ 
rium (Fermi) distribution function for electrons and holes 
in the normal metal and D{pp) is the angle-dependent 
transmission of the tunnel barrier. In the diffusive limit 
the Green functions only weakly depend on the direction 
of the Fermi momentum and, hence, the current (1381) can 
be rewritten in a much simpler form 

I = J - 9ohN + hNgo), (39) 

where gQ’^’^ = indicate the momentum aver¬ 

aged Green-Eilenberger functions and 

^ = e^N^A{\vx\D{pF)) • (40) 

Un 

Combining Eq. (IMl) with the expression for the electric 
potential 

eV = ^ J de{Spg^) = ^ J deSpg^ (41) 

and making use of the condition Sp^,^* = Sp^ff = 0 we 
immediately recover the Ohmic dependence I = V/Rn, 
where the dc photovoltage V should be determined from 
the solution of the quasiclassical equations. Here and 
below the overbar denotes time averaging over the period 
of the external electromagnetic field. 

The task at hand is to evaluate nonequilibrium Green 
functions in a thin superconducting film in the 


E^e^ 

D —^ <C max(T, A), (43) 

this field can be treated as a small perturbation. Within 
this approximation it suffices to disregard the time de¬ 
pendence of the momentum averaged Green functions. 
Replacing Vi i{efc)AiT 3 in Eq. (|A2p and averaging 
this equation over time, we obtain 


X 


Cl{e)+iD —^—^73 {G(£ -I- w) -I- G{s — w)} T3-I- 
{2H(e) - T 3 H(e -|- w)t 3 - T 3 H(e - w)t 3 } , G(e) 


12a;2 
= 0 . 
(44) 


Here we omitted terms proportional to E'^ and for the 
sake of simplicity denoted time averages of the matrices 
G and fl respectively by the same symbols. According to 
Eq. (U) , the matrix H is expressed in terms of the energy 
e, the order parameter A and the self-energy Em- For a 
simple model of randomly distributed isotropic magnetic 
impurities the latter quantity read a^"^P^ 


Sm = +T3U2] ^+igo) ^ 

+ {[ui - T3U2]~^ + igo) (45) 

where riimp is the concentration of magnetic impurities. 
Dimensionless parameters ui and U 2 characterize the 
strength of the spin-isotropic and spin-dependent part 
of the impurity scattering potential. Following our con¬ 
vention we now employ the symbol g^ to denote both 
time- and momentum-averaged Green-Eilenberger func¬ 
tion matrix. 

The relation between the functions g^ and G can be 
derived simply by averaging Eq. (I22p over time. As a 
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result we obtain 


ffo(£) — G{e) + 


2G{e) 


e 

"6 2a;2 

- T3G(e + uj)t 3 - TsGie - oj)f3 


(46) 


Eqs. fully determine a perturbative response of 

a superconductor to an external ac electromagnetic field 

Employing the normalization condition one can param¬ 
eterize the Keldysh component of the Green function ma¬ 
trix in a standard way as 

G^ = G^h - hG^, h = hL+ hhr, (47) 

where the functions and Ht are directly linked to the 
distribution functions for electrons and holes. In equi¬ 
librium we have hL^e) = tanh(e/2T) while the function 
hr equals to zero. In the presence of an ac field the 
function hx already differs from zero and it can be found 
perturbatively from the kinetic equation derived by mul¬ 
tiplying the Keldysh component of the Usadel equation 
by fa and by taking the trace of the resulting ex¬ 
pressions. This equation reads 


Im{A*[E«(s) +F^(5)]} hxie) = -^^^.( 5 ) 

X |i^(£ -I- w) [(i(e) -I- d{e -I- w)] [hLis) - hxie + w)] 

+ i/{e - uj) [d{e) + d{e - Lo)] [hx (e) - /il (e - w)] |, 

(48) 


where d{e) is proportional to the diagonal part of the 
impurity self-energy 


Sp S^(e) - S^(e) = -Aiv{e)d{s). 


(49) 


The function d(e) depends on the impurity scattering 
potential parameters ui ^2 as well as on the quasiparticle 
energy e. We find 

,, s _ u-imp lQuiul{l + ul-ul)lmG^{e) 

“ TtTVo 1(1 + m2 _ ^2)2 + ■ 

Within our model electron scattering on mag¬ 
netic impurities is responsible for charge imbalance 
generationiiii^ as well as for its relaxations^: These two 
effects are accounted for respectively by the left-hand side 
and the right-hand side of Eq. (I48|) . In order resolve this 
equation it suffices to set both the retarded and advanced 
Green functions equal to their bulk equilibrium values. 
The resulting photovoltage value V takes the form 

eV = i y dev{e)hT{e), (51) 



FIG. 2: (Color online) Photovoltage E as a function of the 
microwave radiation frequency ui. The magnetic impurity po¬ 
tential parameters (ui = U 2 = 1.0) and impurity concentra¬ 
tion (riimp = ttNoTco) are the same for all curves. The crit¬ 
ical temperature of a superconductor in the absence of the 
magnetic impurities is denoted as Tco- The inset shows zero 
temperature density of states as a function of energy. 


In Fig. [5] we illustrate typical frequency dependen¬ 
cies of the photovoltage evaluated with the aid of Eqs. 
O-dSlD at different temperatures. The magnetic impu¬ 
rity concentration and the impurity potential parameters 
are chosen to provide the regime of gapless superconduc¬ 
tivity. In this regime the impurity bands overlap with 
the continuous spectrum and the quasiparticle density of 
states remains nonzero in the whole energy range (see 
the inset in Fig. HD. Quasiparticle transitions between 
the states in the continuous spectrum and/or the states 
in the impurity band may yield electron-hole imbalance 
of a different sign depending on the quasiparticle energy 
e. The photovoltage value V is then controlled by an 
interplay between these processes and it may change the 
sign depending on the electromagnetic field frequency. 

In general the voltage E is a complicated function of 
several parameters, such as temperature T, the magnetic 
impurity concentration nimp and the potential parame¬ 
ters Ml and U 2 as well as the frequency of the external 
electromagnetic field uj. Similarly to^ one can expect that 
the maximum value of E is achieved provided both tem¬ 
perature and the magnetic impurity concentration are 
approximately equal to one half of their respective crit¬ 
ical values. In this case and at frequency values several 
times bigger than Tco (cf., e.g., Fig. 2) one can roughly 
estimate 


eE 


Tco 


(52) 


At the border of applicability of Eq. (|i5D one can set 
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E"^ ^ T^Ql{e^VFi) and obtain 

eF ^ T,Ve, (53) 

where Te = tjvp is the elastic electron scattering time. 

V. CONCLUSIONS 

In this paper we extended the well known quasiclassi- 
cal formalism of Usadel equations by deriving the lead¬ 
ing gradient corrections to these equations. At the first 
sight, such corrections should remain small in the diffu¬ 
sive limit and, hence, may at most yield subleading in 
the parameter results. Our analysis demonstrates 
that it is not always so. If, for instance, magnetic im¬ 
purities are present in a diffusive superconductor, spin- 
sensitive electron scattering at such impurities introduces 
electron-hole asymmetry which results in anisotropic in 
the momentum corrections to the quasiparticle distribu¬ 
tion function. Such corrections obviously cannot be cap¬ 
tured within the standard Usadel approximation but they 
are easily recovered with the aid of the formalism devel¬ 
oped here. 

It turns out that the presence of such anisotropic in the 
momentum terms in the distribution function may result 
in dramatic changes in the behavior of a superconductor. 
For instance, the thermoelectric effect in superconduc¬ 
tors (which is usually small in the parameter Tdep) gets 
greatly enhanced in the presence of magnetic impurities 
and the related electron-hole imbalance in the system^. 
We demonstrated that this effect can be conveniently de¬ 
scribed employing the extension of the quasiclassical Us¬ 
adel equation formalism worked out here. 

Our formalism is also well suited in order to describe 
another interesting phenomenon - the so-called photo¬ 


electric effect in superconductors and SN hybrid struc¬ 
tures which implies the appearance of a dc voltage drop 
in the system as a result of inhomogeneous absorption 
of light. Also this effect is strongly enhanced provided 
a certain concentration of magnetic impurities is present 
in a superconductor. If such a superconductor in an SN 
system is irradiated by an external ac electromagnetic 
field the charge imbalance develops which causes a non¬ 
vanishing dc voltage V to occur across a tunnel barrier 
at the SN interface. Here we evaluated the magnitude 
of this voltage and demonstrated that in the presence of 
magnetic impurities in a superconductor the maximum 
value of V is controlled by the electromagnetic signal and 
the electron elastic mean free path i (cf. Eq. (15^ 1 rather 
than by the small parameter Tc/ep which just drops out 
in this case. This large photoelectric effect in SN hy¬ 
brid structures with magnetic impurities can easily be 
detected in modern experiments and can also be directly 
employed in various applications, such as, e.g., a new 
generation of ultrasensitive bolometers. 
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Appendix A 

Making use of Eqs. (fT^ . (ITCl) . (EUl) and (HH) after a 
straightforward but rather lengthy calculation we arrive 
at the following expression for the matrix d which in¬ 
cludes the leading order correction to Eq. m- 


£3 

g^ = -e{GViG - VG - i— (UV,G - HGV, -b GViU - V^GU)-30GVfeGV*GVfeG - 15GVfcGVfcGV*G- 

uf ' 60 I 

-15GV*GVfeGVfcG-llV,GVfeVfeG-llGVfcVfcGV,-b29VfeGVfeV*G-b29GViVfcGVfc-b9GVfeGVfcVi-b9ViVfeGVfcG- 
- llV.GVfcGVfe - llVfcGVfcGV, - 6VfeV,GVfeG - 6GVfcGV,Vfc - 6VfeGV,VfeG - 6GVfcV,GVfc - 6VfeGV,GVfc+ 

-b 6VfeViVfe — GGVfcViVfeG — GV^V^ViG — GViVfcVfcG — V^V^GV^G — GViGVfeVfe -b ViVfcVfc -b 

(Al) 


Employing this result together with Eq. (ED we obtain 
the resulting equation for the matrix function G: 






— — 2ViriVi, Gj -t- i- 


-30V,GVfcGV,GVfe - 29V,VfeGVfcV,- 


Co' VC V t.7j' X80 

+ 15V,GVfcGVfeGV, + 15V,GV,GVfeGVfc + llV,V,GVfeVfc + 6VfcV,GVfcV, + 6VfcGV,VfcV, + 6VfcV,VfcGV,+ 


i'^ k'^ kG'V i + ViGVfcVfeVi + ViGViVfcVfe + V^ViV^GV^, G 

I- 


= 0. (A2) 


The first commutator in Eq. (IA2|) represents the stan¬ 
dard Usadel equation and last two terms define small 
corrections. 

It is easy to observe that within the approximations 
adopted here only the term girii in the expansion (|S]) 
contributes to the electric current 0. Then for the i-th 
component of the current density we have 

ji = - ^^ 24 ^ J dsSp (5/^ ^3) , (A 3 ) 

where is the Keldysh component of the matrix IMI. 


All the equations presented here remain applicable as 
long as the matrix function G changes at distances ex¬ 
ceeding the mean free path i. Provided G changes in 
space sufficiently slowly one can neglect all terms in curly 
brackets in Eq. m as well as the last commutator in 
Eq. (|A2p . Then the above equations receive significant 
simplifications. Further simplifications of our formalism 
occur in the absence of the magnetic field and provided 
the function G depends only one one coordinate, see the 
corresponding discussion in Sec. Ill Cl 
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